Extreme value statistics of networks with inhibitory and excitatory couplings 
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Inspired by the importance of inhibitory and excitatory couplings in the brain, we analyze the 
largest eigenvalue statistics of a random networks incorporating such features. We find that the 
largest real part of eigenvalues of a network, which accounts for the stability of underlying system, 
decreases linearly as a function of inhibitory connection probability up to a particular threshold 
value, after which it exhibits rich behaviors with the distribution manifesting generalized extreme 
value statistics. Fluctuations in the largest eigenvalue remain somewhat robust against an increase 
in system size, but reflect a strong dependence on the number of connections indicating that systems 
having more interactions among its constituents are likely to be more unstable. 
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I. INTRODUCTION 



The largest eigenvalue of network adjacency matrix ap- 
pears in many applications. In particular, for dynamic 
processes on graphs, the inverse of the largest eigenvalue 
characterizes the threshold of phase transition in both 
virus spread |l| and synchronization of coupled oscillators 
[31 in networks. In neuroscience, networks of neurons are 
often studied using models in which interconnections are 
represented by a synaptic matrix with elements drawn 
randomly 0, 3 . Eigenvalues of these matrices are useful 
for studying spontaneous activities and evoked responses 
in such models [3|, |5| , and the existence of spontaneous 
activity depends on whether the real part of any eigen- 
value is large enough to destabilize the silent state in 
a linear analysis. Furthermore, the largest real part of 
the spectra provides strong clues about the nature of 
spontaneous activity in nonlinear models [(|. A recent 
work reveals the importance of the largest eigenvalue in 
determining disease spread in complex networks, where 
epidemic threshold relates with inverse of the largest 
eigenvalue [3]. In context of ecological systems, a cel- 
ebrated work by Robert May demonstrates that largest 
real part of eigenvalue of corresponding adjacency ma- 
trix contains information about stability of underlying 
system ,8]. Mathematically, matrices satisfying a set of 
constraints are stable jOi]. But most real world systems 
have underlying interaction matrix which are too com- 
plicated to satisfy these constraints and hence, study of 
fluctuations in largest real part of eigenvalues is crucial 
to understand stability of a system, as well as of an indi- 
vidual network in that ensemble. 

Largest eigenvalues over ensembles of random Her- 
mitian matrices yielding correlated eigenvalues follow 
Tracy- Widom distribution [l3|. Whereas, extreme value 
statistics for independent identically distributed random 
variables can be formulated entirely in terms of three 



universal types of probability functions: the Frechet, 
Gumbel and WeibuU known as generalized extreme value 
(GEV) statistics depending upon whether the tail of the 
density is respectively power-law, faster than any power- 
law, and bounded or unbounded [ll[ . The GEV statistics 



with location parameter /z, scale parameter a and shape 
parameter ^ has often been used to model unnormalized 
data from a given system. Probability density function 
for extreme value statistics with these parameters is given 

bygU 
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(1) 
Distributions associated with ^ > 0, =0 and < are 
characterized by Frechet, Gumbel, and Weibull distribu- 
tion respectively. Extreme statistics characterizes rare 
events of either unusually high or low intensity. Recent 
years have witnessed a spurt in activities on GEV statis- 
tics, observed in a wide range of systems from quan- 
tum dynamics, stock market crashes, natural disaster 
to galaxy distribution |ll| - [l3| . These distributions have 
been successful in describing the frequency of occurrence 
of extreme events. The experimental examples of GEV 
distributions include power consumption of a turbulent 
flow 1141 ■ roughness of voltage fluctuations in a resistor 
flow |15| . orientation fluctuations in a liquid crystal 16], 
plasma density fluctuations in a tokamak |l7[. Further- 
more, eigenvalues oi a, N x N non-Hermitian random 
matrix with all entries independent, mean zero and vari- 
ance 1/JV, lie uniformly within a unit circle in complex 
plane [18|. Limiting behavior of spectral radius of non- 
Hermitian random matrices has been perceived to lie out- 
side the unit disk as N ^- (X), and with proper scaling and 
shifting, has been found to comply with Gumbel distribu- 
tion [l^ . Though a lot has been discussed about largest 
eigenvalues of random matrices or matrices representing 
properties of above systems, same for adjacency matrices 
of networks has not been done. A vast literature available 
on network spectra is mostly confined to the distribution 
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FIG. 1: (Color online) Largest real part of eigenvalue, av- 
eraged over 2000 realizations of the network, as a function 
of probability of I couplings pi„ for various average degree 
(k) = 10(o), 20(D), 30(0), 40(V), 50(A) and 60(*) from bot- 
tom to top. Left panel is for A'' = 50, middle for A'^ = 100 and 
right panel for A*' — 500. 




of eigenvalues and lower-upper bounds for largest eigen- 
value, etc. [20]. Few available results on the statistics 
of largest real part of network eigenvalues (Rmax) under 
the GEV framework convey that ensemble distribution 
of inverse of Rmax for scale-free networks converges to 
Weibull distribution |21| . Sparse random graphs hav- 
ing N nodes and p connection probability pertains to a 
normal distribution with mean {N — l)p + (1 — p) and 
variance 2p{l — p) [22, [2j]. In the context of brain net- 
works, largest eigenvalues of gain matrices, constructed 
to analyze stability of underlying brain networks, follow 
normal distribution 1241. 



II. RANDOM NETWORK MODEL WITH 
EXCITATORY AND INHIBITORY NODES 

Networks considered in this paper are motivated by 
inhibitory (I) and excitatory (E) couplings in brain net- 
works |25| . entailing matrices with 1,-1 and entries. 
These matrices are different from non-Hermitian random 
matrices studied using random matrix theory framework. 
The entries in matrix for former case take values 0, 1, and 
— 1 instead of Gaussian distributed random numbers. We 
investigate dependence of R^ax on various properties of 
underlying network, particularly on the ratio of I-E cou- 
plings. We find that R-max exhibits a rich behavior as 
underlying network becomes more complex in terms of 
change in I couplings. At a certain I to E ratio, distribu- 
tion manifests a transition to the GEV statistics, which 
is accompanied by another transition from Weibull to 
Frechet distribution as network becomes denser. 

After constructing an Erdos-Renyi random network 
[26] with network size N and connection probability p 
with a corresponding adjacency matrix {A) having entries 
and 1, I connections are introduced with a probability 
Pin as follows. A node is randomly selected as I with 
the probability pm and all connections arising from such 
nodes yield —1 entry into corresponding matrix A. For 
Pin being 0, which assimilates the correlation AijAji — 1, 
network is undirected with A being symmetric entailing 
all real eigenvalues. Maximum eigenvalue for this net- 
work scales as Rmax ~ pN [2^, where quantity pN is 



FIG. 2: (Color online) Statistics of Rmax at different values 
of I coupling probability pi„. The histograms are numerical 
results, solid and dashed lines correspond to the normal and 
the GEV fit respectively. For each case, size of the network is 
N = 100 and connection probability 0.06 which leads to the 
average degree (fc) = 6. (a), (b), (c), (d), (e), (f), (g) and (h) 
correspond to pi„ being 0, 0.1, 0.3, 0.4, 0.42, 0.46, 0.48 and 
0.50 respectively. 



referred to as average degree (fc) of the network. Upon 
introduction of directionality, complex eigenvalues start 
appearing in conjugate pairs, and for pm — 0.5, bulk of 
the eig envalues is distributed in a circular region of ra- 
dius ^yNp{l — p) [28|. Note that for a random network 
with entries 1 and —1, the radius of circular bulk region 
scales with square root of the average degree of the net- 
work i.e. \/pN, and all eigenvalues including the largest 
lie within the bulk. 



III. TRANSITION FROM THE NORMAL TO 
THE GEV STATISTICS 

We investigate Rmax of random networks as a func- 
tion of Pin- Fig. [T] elucidates that, as directionality is 
introduced in terms of I couplings, the mean of Rmax 
decreases linearly up to a certain threshold value, with 
subsequent decrease in a nonlinear fashion without any 
known functional form in terms of network parameters. 
For the linear regime < pi„ < 0.4, fitting with a straight 
line yields the following relation between Rmax and pi„: 



R„ 



pN ~ i2Np)pi, 



(2) 



Fig. [2] depicts statistics for largest real part of eigen- 
value for N = 100 and average degree < k >= 6. The 
curve is fitted with the GEV distribution from Eq. [1] [231 ■ 
For < Pin ^ 0.4, nature of distribution is normal, how- 
ever, as reflected by the left panel of Fig. [TJ the mean 
decreases in agreement to the equation Eq. [5] Variances 
of the data as well as of fitted curves increase with a 
faster rate for < pi„ ^0.1 after which there is a fall 
in its rate of increment. The variance achieves a peak at 
Pin '^^ 0.30, and then decreases with a slower rate. The 
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FIG. 5: (Color onIine)Phase diagram in two parameters space 
A^ and (fc) elucidating nature of GEV statistics based on the 
value of shape parameter ^, and tail of the distribution at 

Pin = 0.5. 



FIG. 3: (Color online) Rmax statistics for network parameters 
N = 100 and p — 0.12 entailing (fc) — 12 average degree, (a), 
(b), (c), (d), (e), (f), (g) and (h) correspond to pi„ = 0.00, 
0.10, 0.30, 0.40, 0.42, 0.46, 0.48 and 0.50 respectively. The 
histograms are numerical results, solid and dashed lines are 
obtained after fitting the data with the normal and the GEV 
distributions respectively. 



behavior of largest eigenvalue statistics is more complex 
in the range 0.40 < pi„ < 0.50, where it can be mod- 
eled using extreme value statistics. Figs, [ll^e)-(h) and 
negative values of the parameter ^ indicate that statis- 
tics converge to the WeibuU distribution. Calculations 
of shape parameter and detailed discussion on fitting has 
been exemplified in the section FVIIII 

As connection probability or (fc) increases, this phe- 
nomena of transition from the normal to the GEV statis- 
tics for Rmax becomes more prominent. Fig. [31 plot- 
ted for N — 100 and < k >= 12, repeats the normal 
distribution behavior for pi„ = 0, which corresponds to 
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FIG. 4: (Color online) For A^ = 100 and p = 0.2 which corre- 
sponds to < fc >= 20 average degree. Subfigures (a), (b), (c), 
(d), (e), (f), (g) and (h) correspond to Pm being 0, 0.1, 0.3, 
0.4, 0.42, 0.46, 0.48 and 0.5 respectively. The histograms rep- 
resent numerical result, solid and dashed lines are obtained by 
fitting the date with the normal and the GEV distributions 
respectively. 



a symmetric random matrix with entries and 1. Till 
Pin ^ 0.3, Rniax Statistics more or less conforms to the 
normal distribution. At pm = 0.4, the statistics deviates 
from the normal distribution, with fitting accuracy being 
higher for the GEV. With a further increase in the value 
of Pin, there is a transition to the GEV statistics as il- 
lustrated by Fig. [3] at pi„ ^ 0.46. This behavior of the 
Rmax continues thereafter. 

As (fc) increases further, p{Rraax) keeps emulating the 
normal distribution at p^ = and the GEV statistics 
at Pin = 0.5. At intermediate Pm values, it manifests 
different behaviors than demonstrated for lower connec- 
tion probabilities as described by Figs. [Hand [31 As soon 
as Pin increases from value 0, the Rmax statistics starts 
deviating from the normal distribution, and for interme- 
diate Pin values, for example at pi„ = 0.2 and pi„ = 0.3 
in Fig. 31 it neither fits with the normal nor with the 
GEV statistics. As value of pi„ increases, Rmax statistics 
indicates a closer fitting with the GEV, and more devia- 
tion from the normal at pin = 0.4 — 0.42 as implied from 
Fig. |4ljd)-(e). Further increase in pi„ prompts a good 
fitting with the GEV statistics at pi„ = 0.44, and this 
good fitting persists thereafter. Detailed discussion on 
true GEV statistics is provided in the section IVIIII 

Aforementioned behavior indicates that smaller (k) 
values induce a smooth transition from the normal to 
the GEV statistics, and for almost all values of pi„ the 
largest eigenvalue statistics remains close to either one of 
them. Whereas larger (fc) values construe a rich behavior 
of p(i?inaa:). It cusucs the normal distribution till certain 
range of pi„ and after that manifests deviation from it 
displaying a rapid change in the statistics as pi„ is incre- 
mented. For this intermediate pi„ range Rmax statistics 
deviates from the normal as well as the GEV substan- 
tially. As Pin increases further, the statistics fits better 
for the GEV as compared to the normal, finally elucidat- 
ing a legitimate fitting with the GEV distribution at pi„ 
being 0.5. 
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FIG. 7: (Color online)Tail behavior of probability density 
function for real part of eigenvalues for network size N — 100 
and Pin = 0.5. Circles represent data points, and solid line 
represents fitting with a straight line. Figures are plotted for 
three different average degrees (a) (fe) = 12, (b) (fc) = 20 and 
(c) (fc> = 50. 



FIG. 6: (Color online) Statistics of Rmax at different values 
of I coupling probability pi„. The histograms are numerical 
results, solid and dashed lines correspond to normal and GEV 
fit respectively. For each case, size of the network is A'' = 100 
and connection probability 0.5 leading to average degree (k) = 
50. (a), (b), (c), (d), (e), (f), (g) and (h) correspond to pi„ 
being 0, 0.1, 0.3, 0.4, 0.42, 0.46, 0.48 and 0.50 respectively. 



IV. TRANSITION FROM WEIBULL TO 
FRECHET 



At these pi„ values where statistic fits well with the 
GEV, the parameter ^, in the tables of section rVlIIl re- 
veals that indeed the distribution complies one of the 
three different statistics, viz. Gumbel, Weibull and 
Frechet, depending upon (k). For small (fc), correspond- 
ing to sparser networks, the GEV statistics espoused 
Weibull distribution, whereas with an increase in con- 
nection probability it indicates a transition to Frechet 
distribution through Gumbel. Phase diagram Fig. [S] il- 
lustrates this behavior for various values of N and (fc). 
For a definite shape parameter range the Weibull and the 
normal states have a close resemblance, the statistics in 
the intermediate regions oipm consequently emulating to 
one of them. Whereas, Gumbel and Frechet are much de- 
viated from the normal, hence in the transition from the 
normal to the Gumbel or Frechet, p{Rmax) niay not abide 
by any of the statistics, and explains a scabrous behavior 
of Rmax in the intermediate pi„ region. For the larger 
(k) values, Rmax does not apprise GEV statistics even 
at Pin — 0.5, Fig. ini and the value of ^ reflect a Frechet 
behavior although KS test rejects it. In order to under- 
stand such an impact of denseness on Rmax behavior, 
we investigate tail behavior of the parent distribution, 
and Fig. [7]^c) reveals that it is deviated from a power- 
law behavior for larger Rmax values, manifesting a devi- 
ation from GEV distribution, whereas tail behaviors cor- 
responding to {k)= 12 and (fc)= 20 imitates exponential 
and power law decay, respectively as indicated by Fig. [7] 
(a) and (b), reinforcing Gumbel and Frechet distribution 
respectively for their maxima. Higher (fc) > 20 values 
for which spectra do not exhibit GEV even for pi„ = 0.5, 
may be ascribed to the correlation in spectra arising from 



1 and —1 entries competing with each other. Fig. [7] indi- 
cates existence of two different scales for (fc) = 50, pro- 
viding a plausible explanation of deviation from GEV. 
Furthermore, revelation of the transition from Weibull 
to Frechet as a function of connection probability or av- 
erage degree (fc) of the network, adds networks to the list 
of wide physical systems exhibiting this transition. For 
example, extreme intensity statistics in relation to com- 
plex random states manifest the Weibull distribution in 
case of minimum intensity^ and the Gumbel distribution 
for maximum intensity [l2j. For mass transport models 
distribution of largest mass displays the Weibull, Gumbel 
and Frechet distribution depending upon critical density. 
|3(]j . For non-interacting Bosons, level density follows 
one of these three distributions depending upon charac- 
teristic ex pon ent of growth of underlying single particle 
spectrum |3l| . 

The interpretation of our result of transition from 
Weibull to Frechet in terms of the stability of underly- 
ing systems can be drawn as follows. For large number 
of I nodes present in the network, the statistics of Rmax 
for denser networks are more right skewed and more de- 
viated from a normal distribution as compared to the 
sparser networks, which indicates that higher values of 
Rmax are more probable for denser networks. This tran- 
spires that the probability with which a network ushers 
to an unstable system is more for denser networks than 
for the sparser ones. Robert May, in his landmark paper 
[8|], concluded that a randomly assembled web becomes 
less robust (measured in terms of its dynamical stability) 
as its connectivity increases. Our results supports this in- 
terpretation for the networks having I and E couplings, 
which is not only based on the average mean behavior 
of largest eigenvalue but also based on its distribution 
modeled using the GEV statistics. 



V. IMPACT OF I-I/E-E AND I-E/E-I 
COUPLINGS 

Our model elucidates a profound impact of I-E ratio on 
both the mean and statistics of Rmax , hence indicating a 
probable impact on the stability or dynamical properties 
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FIG. 8: Plots of Rmax against fintra/finter for diflerent values 
of average degree at pi„ = 0.5 with network size A'" — 100 and 
sample size 20000. Panels (a), (b) and (c) corresponds to k 
= 6, 12 and 20 respectively. 



of corresponding system. To get insight into the transi- 
tion from one statistics to another, first we discuss the im- 
portance of I-E couplings, followed by an analysis based 
on a measure capturing I-E couplings ratio. There exists 
plenty of behaviors and processes exhibited by neural sys- 
tems which have been attributed to the ratio or balance 
between E and I inputs |32l . |33| . In cortex, inter- neurons 
responsible for inhibition play an important function in 
regulating activity of principal cells. When inhibition 
is blocked pharmacologically, cortical activity becomes 
epileptic 3J|, and neurons may lose their selectivity to 
different stimulus (S^- These and other data indicate 
that an interplay between excitation and inhibition por- 
trays a substantial role in determining the cortical com- 
putation [3g. 

In order to understand the origin of two different statis- 
tics at Pin and 0.5, we define a measure fmtra/ finter 
which captures an competition between I-I and E-E cou- 
plings with I-E couplings. The quantities fmter and 
fintra corrcspond to fraction of (I-E) and (I-I) -I- (E- 
E) couplings respectively. Fig. [5] plots Rmax against 
fintra/ finter exhibiting a positivc correlation between the 
two. Presence of few scattered dots towards the right- 
most top corner of the Fig.[S]Ja) for (k) = 6 clearly reveals 
that underlying network has maximum -intra connections 
owing to high fintra/ finter and Rmax- These figures indi- 
cate that the connections between neurons akin escort to 




FIG. 9: Schematic diagram illustrating (a) a modular-type 
structure, (c) a bipartite type structure, and (b) the inter- 
mediate of these two extremes. Dots and arrows represent 
inhibitory and excitatory links, respectively. 



more of an unstable system as compared to a balanced 
structure |6|. Moreover, in realistic neuronal network, 
connectivity is sparser between excitatory neurons than 
between other pairs [37'| , which correspond to the region 
lying towards the left of the Fig. [8] suggesting that net- 
works with less intra-connections are more stable. The 
measure fintra/ finter IS bounded between two extreme 
structures: modular (all I-I and/or E-E connections) and 
bipartite (all E-I or I-E connections) (Fig. |9]). 

Various realizations of the considered model may 
induce networks having (i) modular type structure 
(Fig. [3^), (ii) bipartite type of structure (Fig. [HJ:) and 
(iii) intermediate structure lying in between these two 
(Fig. [3]d). Note that network structure remains same 
in all three cases, it is only the type of node (I or E) 
at two ends of a connection which decides the config- 
urations mentioned above. An ideal bipartite structure 
would bring upon an anti-symmetric matrix consequently 
having all imaginary eigenvalues. Though networks con- 
sidered here do not consort to an ideal bipartite arrange- 
ment as depicted in Fig. [SJc), for high values of finter as 
elucidated in the Fig. [H it is expected to lie close to this 
arrangement which explains the origin of lower R-max to 
the left of Fig. |51 What follows is that larger I-E cou- 
plings drives to lower values of Rmax , which may be even 
for an ideal case of bipartite structural arrangement en- 
tailing a complete anti-symmetric matrix, whereas larger 
I-I or E-E couplings, which may be considered as modu- 
lar type arrangement direct to higher Rmax values which 
may sometimes be unusually very large for certain net- 
work configurations, probably being one of the plausible 
reasons behind the origin of GEV statistics. Further- 
more, the discussion elaborating Fig. |H] apparently sheds 
light on the origin of stability of network configurations 
having more inter-connections, in turn supporting bipar- 
tite type topology over a modular one as proposed in [3^ 
for real world network. 



VI. CONCLUSION AND IMPLICATIONS 

To conclude, we have analyzed Rmax statistics of net- 
works having I and E couplings. A linear decrease, fol- 
lowed by a non-linear one, in R-max as a function of pin 
indicates that an increase in complexity, in terms of in- 
clusion of I nodes, increases the stability of underlying 
system. For the range where R-max niean ensues the lin- 
ear dependence on pi„, the statistics mostly yields a nor- 
mal distribution, and after this critical pi„ value there 
is a transition to the GEV statistics. The versatile sit- 
uation arising from I-I, I-E competition, bringing upon 
GEV statistics, has not been observed for zero Pm value, 
and hence may be attributed to the rich behavior of R-max 
in the presence of I nodes. 

Though modeling real brain networks needs to account 
for more properties such as specific degree distribution, 
hierarchical structure etc, which may bring upon a richer 
largest eigenvalue pattern [39| , an impact of I nodes im- 



pels a drastic change in its spectral properties illustrating 
extreme events which has been envisaged upon in this pa- 
per. Asymmetric matrices considered here, motivated by 
brain networks, elucidate a different statistical property 
of Rraax than that of non-Herniitian matrices motivated 
by ecological webs [401 ■ Moreover, the universal GEV 
distribution displayed by largest eigenvalue of networks 
propagates theory of extreme value statistics, which sug- 
gests that a model which fits with all eigenvalues or de- 
scribe fluctuations of all eigenvalues 4l| may not be a 
good model for the largest one. 

Recent years have seen a fast development in merg- 
ing of extreme statistics tools and random matrix the- 
ory. The present work extends this general perspective 
to complex networks. To our knowledge, this is the first 
work on networks demonstrating that the largest eigen- 
value of a network, at particular I-E coupling ratio, can 
be modeled by the GEV statistics. The transition of 
the statistics from one type to another as a function of 
I connections has crucial implications in predicting and 
analyzing network functions and behaviors in extreme 
situations [42.] . 



VII. ACKNOWLEDGMENT 

SKD acknowledges UGC for financial support. SJ 
thanks DST for funding. It is a pleasure to acknowledge 
Dr. Changsong Zhou (Hongkong Baptist University) for 
useful discussions on brain networks at several occasions, 
and Dr. A. Lakshminarayan (IITM) for suggestions on 
GEV statistics. 



VIII. APPENDIX 

We use Kolmogorov-Smirnov (KS) test to characterize 
hypothesized model of our data. The KS test is known 
to be superior to other techniques such as chi-square test 
[43| for identifying a particular distribution. For exam- 
ple, in the context of networks, the said test has been 
performed to confirm power law for a given network data 
IHl. The function kstest of MATLAB Statistics Toolbox 
is used to verify the acceptance of a given statistics at 
95% level of confidence if its corresponding p- value of KS 
test is greater than 0.05. 

In some of the parameter regimes, GEV distribution 
resembles the normal distribution owing to its shape pa- 
rameter ^, which characterizes it as Weibull distribu- 
tion [45i] ■ and a particular distribution is confirmed us- 
ing KS test. Another example demonstrating the quality 
of our results can be exemplified with larger (fc) values, 
where for pm > 0.46, though distribution looks more like 
Frechet (Fig.4(d)), KS test accepts Frechet distribution 
at Pin — 0.5 only. We perform KS test for sample size 
5000, which is large enough to approve a statistics. For 
example, 46] accounts for 4000 sample size for perform- 
ing KS test, and in [43], it is implemented to affirm GOE 



TABLE I: Estimated parameters of GEV and normal distri- 
butions for Rmax for different inhibitory coupling probability 
(pin). For each case, size of network is A'^ = 100 and average 
degree (fc) = 6. Sample size is 5000 for all pi„ values except 
for • entries for which sample size is 20000. 



Pin 


e of 
GEV 


a of 
GEV 


/i of 
GEV 


P- 

value 

of KS 

test 

for 

GEV 


fl of O" of p- 
Normal Normal value 
of KS 
test 
for 
Norma' 


0.00* 


-0.2392 


0.3509 


6.8484 


0.0013 


6.9813 


0.3548 


0.3717 


0.10 


-0.2204 


0.4806 


5.8032 


0.0003 


5.9893 


0.4800 


0.4127 


0.30* 


-0.2248 


0.5410 


4.0107 





4.2210 


0.5505 


0.2966 


0.40 


-0.1945 


0.5230 


3.2444 


0.0062 


3.4606 


0.5501 


0.0000 


0.42* 


-0.1695 


0.4960 


3.0695 


0.0637 


3.2852 


0.5355 


0.0000 


0.46 


-0.0933 


0.4178 


2.8485 


0.3270 


3.0558 


0.4832 


0.0000 


0.48 


-0.0881 


0.3767 


2.7845 


0.3983 


2.9725 


0.4374 


0.0000 


0.50 


-0.1104 


0.3492 


2.7593 


0.9919 


2.9261 


0.3955 


0.0000 



TABLE II: Estimated parameters of GEV and normal distri- 
butions for Rmax for different inhibitory coupling probability 
(pin). For each case, size of network is N = 100 and average 
degree (fc) = 12. 



Pin 


e of 
GEV 


a of 
GEV 


/i of 
GEV 


P- 

value 

of KS 

test 

for 

GEV 


H of cr of p- 
Normal Normal value 
of KS 
test 
for 
Norma 


0.00 


-0.2482 


0.4693 


12.611 


0.0262 


12.7853 


0.4702 


0.6673 


0.10* 


-0.3009 


0.8671 


10.272 





10.566C 


0.8475 


0.0001 


0.30* 


-0.2370 


1.0494 


6.1683 


0.0187 


6.5675 


1.0617 


0.4473 


0.40 


-0.1742 


0.9214 


4.4642 


0.0136 


4.8629 


0.9946 


0.0001 


0.42 


-0.1135 


0.8319 


4.1918 





4.5937 


0.9456 


0.0000 


0.46 


0.0584 


0.5809 


3.7219 


0.0001 


4.0940 


0.7793 


0.0000 


0.48 


0.0711 


0.4908 


3.5951 


0.0709 


3.9159 


0.6771 


0.0000 


0.50 


0.0232 


0.4500 


3.5744 


0.7076 


3.8454 


0.5916 


0.0000 



and GSE statistics for random matrices with sample size 
1000. 

It might be possible that for some network parameters, 
KS test accepts normal as well as Weibull distributions, 
as depicted earlier by the fact that GEV distribution in 
a certain shape parameter range resembles normal distri- 
bution '49] . To address this issue, we increase the sample 
size from 5000 to 20000 for which KS test accepts either 
normal or Weibull distribution. For example at (fc) = 6 
for various pm values 0.0, 0.3, 0.4, the sample size is in- 
creased to 20000 where only one distribution is accepted 
by the KS test. Similarly for (fc) = 12 and pi„ = 0.1 and 
0.3, the sample size is increased to 20000 for implemen- 
tation of KS test. For (fc) values ranging between 16 
and 20, the distribution lies close to Frechet but not ex- 
actly Frechet even for 20000 sample size, thus rendering 



TABLE III: Estimated parameters of GEV and normal distri- 
butions for Rmax for different inhibitory coupling probability 
(pin)- For each case, size of network is N = 100 and average 
degree (fc) = 20. 



Pin 


C of 
GEV 


a of 
GEV 


H of 
GEV 


P- 

value 
of KS 

test 

for 

GEV 


fl of (T of p- 

Norma] Normal value 
of KS 
test 
for 
Normal 


0.00 


-0.2293 


0.5546 


20.388 


0.0443 


20.601 


0.5587 


0.9446 


0.10 


-0.2999 


1.3779 


16.325 


0.0025 


16.790 


1.3371 


0.0183 


0.30 


-0.2420 


1.8081 


8.7824 


0.0037 


9.4622 


1.8091 


0.2053 


0.40 


-0.1231 


1.4617 


5.7608 





6.4558 


1.6475 


0.0000 


0.42 


-0.0221 


1.2347 


5.2248 





5.9203 


1.5168 


0.0000 


0.46 


0.1984 


0.7740 


4.5096 





5.1254 


1.2005 


0.0000 


0.48 


0.1836 


0.6413 


4.3695 


0.0035 


4.8724 


1.0147 


0.0000 


0.50 


0.1133 


0.5494 


4.3218 


0.1123 


4.7080 


0.8266 


0.0000 



KS test to reject it. This is supposedly the bottleneck 
of increasing sample size. We perform KS test for even 
a higher sample size (50000), and it does not accept the 
Frechet distribution (even though distribution keeps ly- 
ing close to the Frechet), hence demonstrating fairness of 
our data and the technique adopted to conclude a par- 
ticular distribution. 

We have also observed an effect of network size on the 
value of shape parameter ^. For example (fc) = 6, net- 
works size TV = 100 and N = 1000 yield a ^ which char- 
acterizes WeibuU distribution, whereas for (fc) = 20, size 
N = 100 reflects Frechet, and size N = 1000 reflects 
Gumbel distribution. The phase diagram presented in 
Fig. [5] corresponds to pi„ = 0.5, for which we get GEV 
statistics till certain (fc) values. For the larger (fc) values 
when Rinax does not comply with GEV statistics even at 
Pin =0.5, Fig. [6] and the value of ^ in the Table. 4 suggest 
a Frechet behavior however KS test rejects it. 



TABLE IV: Estimated parameters of GEV and normal distri- 
butions for Rmax for different inhibitory coupling probability 
(pin). For each case, size of network is N = 100 and average 
degree (k) = 50. Since for none of the pi„ values data fits 
with the GEV distribution, sample size here is increased from 
5000 to 20000 for all pi„ values to inquire if higher sample size 
confirms a GEV distribution. 



Pin 


i of 

GEV 


a of 
GEV 


/i of 
GEV 


P- 

value 
of KS 

test 

for 

GEV 


/x of 
Norma 


a of 
Normal 


P- 

value 

of KS 

test 

for 

Normal 


0.00 


-0.2354 


0.7033 


49.714 


0.0000 


49.982 


0.7073 


0.9484 


0.10 


-0.2787 


3.2684 


38.806 


0.0000 


39.955 


3.1831 


0.0000 


0.20 


-0.2803 


4.2905 


28.419 


0.0000 


29.920 


4.1866 


0.0004 


0.30 


-0.2304 


5.0941 


17.956 


0.0000 


19.897 


5.0239 


0.0000 


0.40 


0.0676 


3.4118 


8.2755 


0.0000 


10.489 


4.4492 


0.0000 


0.42 


0.4561 


2.2051 


6.6802 


0.0000 


9.0732 


3.9636 


0.0000 


0.44 


0.5707 


1.5312 


5.9710 


0.0000 


7.9807 


3.3925 


0.0000 


0.46 


0.5063 


1.1223 


5.6447 


0.0000 


7.1094 


2.7603 


0.0000 


0.48 


0.3859 


0.8741 


5.4775 


0.0000 


6.4787 


2.1008 


0.0000 


0.50 


0.24152 


0.7562 


5.4673 


0.0000 


6.1574 


1.5986 


0.0000 
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